Competing interactions in arrested states of colloidal clays 
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Using experiments, theory and simulations, we show that the arrested state observed in a colloidal 
clay at intermediate concentrations is stabilized by the screened Coulomb repulsion ( Wigner glass) . 
Dilution experiments allow us to distinguish this high-concentration disconnected state, which melts 
upon addition of water, from a low-concentration gel state, which does not melt. Theoretical 
modelling and simulations reproduce the measured Small Angle X-Ray Scattering static structure 
factors and confirm the long-range electrostatic nature of the arrested structure. These findings are 
attributed to the different timescales controlling the competing attractive and repulsive interactions. 
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Dynamical arrest in soft colloidal systems has recently 
become the subject of an intense research activity. The 
fine tuning of control parameters opens the possibility 
to tailor the macroscopic properties of the resulting non- 
ergodic states. Several mechanisms of dynamical arrest 
have been identified. Building on the knowledge on the 
hard spheres glass [ij, it has become recently clear that 
when both attractive and repulsive terms are present in 
the interaction potential, a re-entrant liquid-glass line, 
surrounded by two distinct glasses, has been predicted 
and experimentally observed in short-ranged attractive 
colloids at high concentrations 2]. A rich phenomenol- 
ogy also takes place at low concentrations: here gelation 
occurs, which may result from different routes [3|. Inter- 
esting scenarios arise when, in addition to a short-ranged 
attraction, particles have a residual electrostatic charge 
which builds up a long-range repulsion in the effective po- 
tential. In this case, particles can form equilibrium clus- 
ters [3], which provide the building blocks of arrest 3. 
Recent works have shown that both Wigner glasses 6|, 
intended as arrested states formed by disconnected par- 
ticles or clusters and stabilized by the electrostatic repul- 
sion, and equilibrium gels, which occur at larger packing 
fractions when the clusters branch into apercolating net- 
work, can form under these conditions 



To investigate the formation of multiple arrested 
states, colloidal clays have emerged as suitable 

candidates. The anisotropy of the particles, combined 
with the presence of attractive and repulsive terms in 
the interactions, makes the phase diagram of such col- 
loidal systems very complex. Among these, Laponite 
suspensions have been widely studied not only for their 
appealing industrial applications [lT| but also for their in- 
teresting experimental/theoretical properties [l2l - Eoj |. In 
particular, Laponite displays a non-trivial aging dynam- 
ics [iBl and (at least) two final arrested states, which are 
obtained by a simple increase of Laponite volume frac- 
tion from low {Cw < 2.0%) to moderate {C^ > 2.0%) 
values, at fixed salt concentration Cg = 10~^M~[l5, 17 1. 



More recently, the static properties of these two states 
have been investigated in detail 18], showing that they 
are also characterized by both the very different aging 
behaviors and the shape of the static structure factor 
S{Q). While the low C^, state shows typical inhomo- 
geneities, the high Cw state is homogenous, allowing for 
the identification of these two arrested states respectively 
as gel and glass. 

In this Letter, through the combination of experi- 
ments, numerical and theoretical approaches, we demon- 
strate that the non-ergodic state observed in Laponite 
samples in salt free water at sufficiently large concentra- 
tions {Cw > 2.0%) is a Wigner glass. By performing 
a simple but impressive dilution experiment, we are able 
to distinguish whether attractive or repulsive interactions 
are dominant in the formation (and stability) of the ar- 
rested structure for both low and high concentration ar- 
rested samples. Moreover, we compare the S{Q) deter- 
mined by Small Angle X-rays Scattering (SAXS) with 
theoretical and numerical ones in the high concentration 
window, providing an estimate for the effective interac- 
tions between Laponite platelets. These results provide 
clearcut evidence that the resulting large concentration 
arrested state (which is still found at very low packing 
fraction) can be attributed to the residual electrostatic 
repulsion (Wigner glass) 6] . Thereby Laponite offers the 
unique example of a system which, not only displays two 
distinct non-ergodic states differing only in colloidal con- 
centration, but also shows a counter-intuitive scenario 
which involves gel-like structures at lower concentrations 
and a (truly) disconnected glass at larger ones. Hence, 
the increase of concentration, which normally favours 
percolation and branching, here has the opposite effect. 
We attribute this novel behavior to the interplay of dif- 
ferent time-scales in the involved attractive and repulsive 
interactions, a feature that may be found in other com- 
plex colloidal systems. 

To investigate the glassy state formed by Laponite at 
high clay concentration (C„ > 2.0%) [ll, [H [3, we 
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FIG. 1: Evolution of a high concentration sample at Cw = 3.0% with time. The initially liquid (L) sample (A) ages with time 
reaching, after 50 h, a final arrested glassy state (G2)(B). At this time a colored solution (in a liquid (L) state) is added to the 
arrested sample (C) that starts to fluidize at the interface (D)- 70 h. Finally the whole sample is liquid again (E) - 120 h. 



use a combination of different approaches. Firstly, we 
perform a dilution experiment (adding water to the ar- 
rested state) that provides a direct criterion to establish 
the dominant interactions controUing the formation and 
stability of the arrested state 0]. If this were due to 
attractive interactions, the presence of additional water 
should not affect the final state, because Laponite bonds 
(of strong electrostatic nature) can not be broken by the 
mediation of water. On the other hand, if repulsive inter- 
actions are determining the formation of the non-ergodic 
state, the increase of free volume allows a rearrangement 
of the Laponite platelets on average to a larger distance. 
If such distance is larger than the characteristic repulsive 
length (namely the electrostatic Debye length) the glass 
will be destabilized up to melting to a liquid state. 

To obtain reliable and reproducible results 21 1 
Laponite RD dispersions are prepared using the same 
protocol described in [18[. The starting waiting time 



[t^ — 0) is the time at which the suspension is filtered. 
For the dilution experiment a sample at high clay concen- 
tration, Cyj = 3.0%, is filtered and sealed in small glass 
bottles. The photographic sequence of this experiment 
is shown in Fig[TJ The initially liquid (L)((A) panel) 
sample ages with time « 50 h [15| until reaching a final 
arrested state, as it can be seen in panel (B) where the 
sample shows its solid like nature (G2). At this point 
a solution of colored water (of the same volume as that 
of the Laponite solution) is added to the arrested (col- 
orless) sample. The colorant used, Rhodamine B at a 
concentration 10"'^ Af, does not play any role in the pro- 
cess, as shown below, and is used to better distinguish (C 
panel) the arrested sample (G2) from the liquid solution 
(L). After the addition of water, the arrested sample 
starts to fluidize at the interface with the liquid solution 
(D panel - — 7Q h). This process evolves with time 
until, after 120 h, the whole sample is back into a liq- 
uid state (E panel). This result strongly speaks in favor 
of the repulsive nature of the Cw — 3.0% arrested state. 
The experiment is repeated in absence of Rhodamine and 
shows that the addition of the dye does not affect the re- 
sults. We have monitored that the pH of the solution 
remains constant throughout the dilution experiment. 



We now turn to elucidate the effective interactions be- 
tween Laponite platelets, by comparing the evolution of 
the experimental static structure factor 5*^((3), obtained 
from SAXS measurements performed at the High Bril- 
liance beam line (ID2) [22i] at the European Synchrotron 
Radiation Facility (ESRF) in Grenoble, France, as de- 
scribed in [1^, with that predicted theoretically S*^{Q), 
assuming several effective interaction potentials [2^ . 

In Fig. [2] the comparison between 5*^((5) (symbols) 
and optimal S^^{Q) (lines), for fixed t^, ~ 50 h in the 
concentration window 2.0% < Cw < 3.0%, is shown. As 
it is evident from Fig. [21 as concentration is decreased 
from Cw = 3.0% to = 2.0% the low Q signal is slightly 
increased and the main peak is slightly shifted towards 
lower Q values. The position of the main peak is around 
Q ~ 0.17 nm^^, corresponding to a length w 37 nm, a 
value clearly larger than the platelet diameter d = 25 nm, 
pointing to a non-connected final structure. Moreover, no 
evidence of a peak at contact distance is found. 

The theoretical S^'^{Q) is calculated by solving numer- 
ically the Ornstein-Zernike equation with an appropri- 
ate closure relation or by direct simulation of the 
effective potential, without a significant change in the 
description of the experimental data. For simplicity, we 
report Percus-Yevick results in the following. In this sim- 
plified treatment, we consider only the centers of mass of 
the scattering objects, supposedly single platelets, and 
treat them as (spherical) points, following previous work 
on different systems Q- The main fit parameters in- 
volved in the description of the experimental data are 
the number density p of the scattering objects and the 
number of counter-ions Nd that dissociate (on average) 
from each platelet when dissolved in water. Our study 
indicates that p cannot be the same as the number den- 
sity of Laponite discs obtained from the nominal weight 
concentration p^ — 2.4 • 10^^/m^. Indeed, no simple in- 
teraction potential is able to reproduce the position of 
the main peak using p^. Hence p must be smaller, due 
to the fact that plat elets may be found within a distribu- 
tion of clusters [25[. We find that the data can be fully 
described, in terms of peak position, peak height and 
compressibility, by a Yukawa potential, which accounts 
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FIG. 2: Comparison between measured S^' (Q) (symbols) 
and theoretical S^'^{Q) (lines) for several high concentration 
Laponite samples at ~ 50 h. For clarity, curves have been 
shifted along the vertical axis progressively by 0.3. Inset: 
Comparison between S^'^ (Q) (symbols) and S^'^{Q) (lines) for 
two different samples with Cw ~ 3.0%. 



in an average way [tI for the screened electrostatic repul- 
sion between Laponite platelets or clusters, in a gree ment 
with more sophisticated theoretical approaches [19|. We 
use a potential of the form Fg// = A^e~^^^ /r, where ^ is 
the Debye screening length and A provides a measure of 
the strength of the repulsion. The number density is fixed 
to p ~ 0.38pi, a relationship which is found to hold in the 
whole investigated concentration window. This suggests 
that Laponite may be dispersed in a distribution of very 
small aggregates, possibly due to particles that are not 
completely delaminated, compatibly with an AFM study 
performed under very dilute conditions [2^. The num- 
ber of counterions in solution is optimized to Nd = 60, 
a value well below the bare charge of a single platelet 
(700 e), in good agreement both with previous sim- 
ulations [19I] and with conductivity measurements [TtJ. 
With this choice of parameters, the change in nominal 
concentration and its associated screening length ^, vary- 
ing between 8 and 10 nm, are able to reproduce S^^{Q) 
in the full investigated window, as seen in Fig. [2j The 
repulsion strength A^/ksT is found to increase in the 
studied range of concentration by a factor of « 2, com- 
patibly with the behavior of other charged systems [i^ . 

While the Yukawa potential correctly reproduces the 
experimental data, other possible candidates do not. The 
simple hard sphere model does not capture at all the 
softness of the interactions, failing to reproduce the low- 
Q regime also after one adjusts p. On the other hand, 
the addition of a short-range attraction to the electro- 
static repulsion is (i) not sufficient to explain the low 
Q-position of the main peak; (ii) tends to induce an en- 
hanced structuring at nearest-neighbour length. There- 
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FIG. 3: Comparison between simulated (symbols) and theo- 
retical (lines) S(Q) for Cw = 2.4% and 3.0% Laponite sam- 
ples. Inset: Simulated (cross symbols), theoretical (line) and 
measured (open circles) S{Q) for Cm = 2.4%. 



fore, we conclude that no relevant attractive features are 
present in S{Q) in this concentration window. Hence, the 
combination of the dilution experiment with the theoret- 
ical analysis allows to identify the observed arrested state 
as a Wigner glass stabilized by the residual electrostatic 
repulsions. 

One further comment concerns the robustness and re- 
producibility of our findings for different samples. Al- 
though different measurements can bring small changes 
in S'^'^Q) due to differences in Laponite batches, to the 
filtration procedure or to real clay concentrations, we can 
fit equally well different samples with the same effective 
Yukawa potential upon variation of the fit parameters. 
This is shown in the inset of Fig. [2] where S^^ {Q) is re- 
ported for two different samples with — 3.0%, from 
different batches, measured in two different runs, in com- 
parison with S*'^{Q). For the second batch, which has a 
higher polydispersity, we find Nd = 90 and p = 0.26pL- 

The theoretical approach we have adopted so far treats 
only the correlations between the centers-of-mass of the 
scattering objects. Although such an approach was 
proven to be a good (although crude) approximation to 
describe the static correlations between equilibrium clus- 
ters of non-spherical shape 0, it is crucial to show that 
this treatment is valid also when we take into account 
explicitly the actual disc shape of clay particles. To this 
aim, we have performed MonteCarlo simulations of 200 
discs at the same p as that extracted from the fits. We use 
the interaction model introduced in [23| where the total 
negative charge is uniformly distributed over the surface 
area, while rim charges are neglected. Each Laponite 
platelet is schematized as a rigid disk, composed by 19 
sites disposed on a regular mesh. We simulate several 
repulsion parameters and find that the numerical S"'{Q) 
is in good agreement with the experimental and theoret- 
ical ones when the total charge is fixed to 70 e and the 
screening length is 5 nm. Results are shown in Figure 
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FIG. 4: Starting (A), intermediate (B) and final (C) states of 
a low concentration sample at Cw = 1.5%: the initial gel state 
Gi is not macroscopically affected by the addition of water 
(C). Starting (D), intermediate (E) and final (F) states of a 
high concentration sample at Cw ~ 3.0%: after the addition 
of water (L) the original glassy state G2 fluidizes (F). 



|3] for two samples with = 2.4% and Cw = 3.0%, in 
comparison with S*^{Q). In the inset, simulated, the- 
oretical and measured S{Q) for sample Cw = 2.4% are 
compared. Hence, the favorable comparison of simplified 
theory and simulations with experimental measurements 
provides evidence that our effective approach is useful to 

fite^^l]?t&iil^Srci!l§Bo^ 1 _ 1 o\Kg4. .. 1/22/: 

So far we have reported evidence of the existence of a 

Wigner glass state in Laponite by the combination of ex- 
periments, theory and simulation for 2.0% < C^ < 3.0%. 
At lower concentrations structural measurements have 
shown a marked growth of S{Q) at low Q with waiting 
time, which was attributed to the presence of attractive 
interactions determining an arrested state of attractive 
(or gel) nature [ill ■ Already for the lowest concentration 
shown here {Cw = 2.0%), the agreement between S^^ (Q) 
and S^^{Q) is less good in the low Q range, as visible 
in Fig. [21 This might be an indication of the increasing 
role of attractive interactions upon lowering concentra- 
tion. Hence, we perform the dilution experiment also for 
a low concentration sample expecting a different result 
with respect to the high concentration one. 

In Fig. 2] the dilution experiment is shown for a low 
{Cw = 1.5%) and a high {Cw = 3.0%) concentration 
sample. The samples are left to age up to the final cor- 
responding arrested states, respectively of gel (Gi) and 
glass (G2) nature for low and high concentrations [l8| . 
While the sample at Cw — 3.0% arrests within 50 h, 
1.5% takes several weeks 



the one for Cw = 1.5% takes several weeks [15|, |17[. Af- 
ter arrest takes place deionized water (L) is added and 
panels (A) and (D) show the identical situation for the 
two cases. However the evolution of the two samples 
soon becomes dramatically different. The arrested state 
of the high concentration sample (G2) starts to fluidify 
((E) panel) up to reach a final liquid state ((F) panel). 
On the contrary, the low concentration sample does not 



liquify and does not show any macroscopic changes of 
its solid like state, even after waiting additional weeks, 
i.e. the same timescale of the arrest process. In this 
case the presence of available free volume does not de- 
termine any change in the gel-arrested state ((B) and 
(C) panels), due to the fact that water cannot break clay 
bonds. The different behavior between the high and low 
concentration samples ((C) and (F) panels) observed in 
the dilution experiment confirms the different nature of 
the corresponding arrested states: a Wigner glass (G2) 
and a gel (Gi), dominated respectively by repulsive and 
attractive interactions at high and low concentrations. 

The striking result of a low-concentration gel network 
and of a high-concentration disconnected Wigner glass is 
at odds with previous studies of simpler systems j7i.l8l.l28|. 
We attribute this unexpected scenario to the separate 
timescales controlling the interactions and the two arrest 
processes [lB|. While repulsion is felt almost immedi- 
ately after samples are prepared, attraction, due to its 
anisotropic nature and to the presence of an effective 
repulsive barrier, develops on a much longer timescale 
as in reaction-limited aggregation. Long-time attraction 
may also affect the repulsive Wigner glass, through the 
formation of subsequent additional bonds. Indeed, rhe- 
ological measurements for rejuvenated samples show a 
dependence on the idle time waited after the glass is 
formed [2^. To elucidate this point we have performed 
an additional dilution experiment on the high concentra- 
tion sample waiting one week after dynamical arrest takes 
place. In this case the system does not melt, probably 
due to intervening long-time attraction which strength- 
ens the glass. However the unmelted final state is dif- 
ferent from that observed at low concentration, while 
the latter remains unchanged and water does not pen- 
etrate into it, the former seems to slowly absorb water in 
the bulk. SAXS measurements show that S{Q) does not 
change significantly in this time window demonstrating 
that the disconnected structure of the repulsive glass is 
preserved despite the intervening attraction, and hence 
the resulting state cannot be considered an attractive 
glass in the most common sense. 

In conclusion, through the combination of dilution ex- 
periments, SAXS experiments, theory and simulations, 
we have shown the presence of a disconnected Wigner 
glass state in a screened charged colloidal system. This 
arises at a larger (but still very low) concentration with 
respect to the one where a gel network is found, thanks to 
the different timescales controlling the repulsive (short- 
time) and attractive (long-time) interactions. We expect 
these findings to be relevant to other complex liquids with 
competing interactions as anisotropic patchy colloids [30l | 
and globular proteins [3]| . 
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